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Abstract 

We investigate the production of the novel P-wave mesons D* and D[ (D x ), identified as J p = 0+ 
and 1 + , in heavy B meson decays, respectively. With the heavy quark limit, we give our modelling 
wave functions for the scalar meson Dq. Based on the assumptions of color transparency and 
factorization theorem, we estimate the branching ratios of B — * Dq-tt decays in terms of the 
obtained wave functions. Some remarks on productions are also presented. 
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I. INTRODUCTION 



It is doubtless that quark model provides a successful method to describe the hadron 
spectroscopy. For instance, based on the SU(3) flavor symmetry, quark and anti-quark can 
comprise octet and singlet states, called nonet together, such as the well known pesudoscalar 
mesons of pion, kaon and eta with J PC = CT + . However, if we apply the same concept on 
the light scalar mesons described by ++ , such as the nonet composed of isoscalars cr(600) 
and /o(980), isovector ao(980) and isodoublet k, there are some puzzles (a) why do (980) and 
/o(980) are degenerate in masses and (b) why the widths of o and k are broader than those 
of a (980) and / n (980) [lj. It is probable that these scalar states consist of four-quark rather 
than two-quark |2| . Moreover, the possibilities of KK molecular states, gluonium states and 
scalar glueballs are also proposed. It is clear that the conclusion is still uncertain. 

Now, the mysterious event happens not only in the light scalar meson system, but also in 
the heavy cs one. Recently, BABAR collaboration has observed one narrow state, denoted 
by D*j(2317), from a Dfn° mass distribution j^j]. Furthermore, the same state has been 
confirmed by CLEO and a new state -D*j(2463) is also observed in the D* + tt° final state j^J. 
Finally, BELLE verifies the observations [5|. By the data analysis, D*j(2317) and D*j(2463) 
are identified as parity-even states with + and 1 + , respectively. From the observations, 
the interesting prob.cm is that the states of ^(2317) and ^(2463) cannot match with 
theoretical predictions [6(, i.e., the masses (widths) are too low (narrow). To explain the 
discrepancy, either the theoretical models have to be modified |ZJ or the observed states 
are the new composed states. To satisfy the latter, many interesting solutions have been 

suited rcent™. flayyyyQ. 

In fact, before the BABAR's observation, CLEO [15[ and BELLE [16| measured the 
similar states in the cq (q = u, d) system in heavy B meson decays. With two-quark picture, 
there are four parity-even (angular momentum I — 1) states described by J p = + , 1 + , 
1 + and 2 + , respectively. J = j q + Sq is the total angular momentum of the corresponding 
meson and consists of the angular momentum of the light quark, j q , and the spin of the 
heavy quark, Sq, where j q — S q + £ is combined by the spin and orbital angular momenta 
of the light quark. In the literature, they are usually labeled by D$, D[, D\ and D\, 
respectively. We also use D** to denote all of them. The first two belong to j q = 1/2, 

0(1) ^1(2) 



while the last two j q = 3/2. In the heavy quark limit, it is known that D Q )A and D 
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decay only via S and D-w&ve, respectively. Therefore, one expects that the widths of the 
former are much broader than those of the latter, which is consistent with the observations 
of CLEO and BELLE Q- 

Even the BELLE's updated data [17( also show the same 
phenomenon. We now summarize the results of CLEO and BELLE as follows: in CLEO, 
the masses (widths) of P-wave states are given by m^^r^J = 2422.0 ± 2.1 (I8.9I35) an< ^ 
nT'D'O^D*) = 2458.9 ± 2.0 (23 ± 5) MeV, and the measured branching ratios (BRs) of B 
decays are given as 

BR(B- -> D°n-) x BR(D° -> D* + n-) = (7.8 ± 1.9) x 10" 4 , 
BR(B- -> D*°7T') x BR(Df -> D* + tt") = (4.2 ± 1.7) x 10~ 4 . (1) 

In BELLE, the four states are all measured as m D *(T D *) = 2308 ± 17 ± 15 ± 28 (276 ±21 ± 
18 ±60), m D{ (T D{ ) = 2427.0 ±26 ±20 ±15 (384+^> 7 ± 24 ± 70), m Dl (T Dl ) = 2421.4±1.5± 
0.4±0.8 (23.7±2.7±0.2±4.0), andm D *(r D *) = 2461.6±2.1±0.5±3.3 (45.6±4.4±6.5±1.6) 
MeV, and the BRs with possible decaying chain being 



BR(B- -> D*°n-) x BR(D*° -> P+tT 
BR(B~ -> DjV - ) x BR(D? -> P*+tT 

PP(P~ -> P°7r~) X Ml/)? -> P* + tT 



PP(P~ -> D^tt") x BR(Df -> P 



7T 



PP(P~ -> P^O x BR(D? -> P* + 



7T 



(6.1 ±0.6 ±0.9 ±1.6) x 10~ 4 , 
(5.0 ±0.4 ±1.0 ±0.4) x 10" 4 , 
(6.8 ± 0.7 ± 1.3 ±0.3) x 10~ 4 , 
(3.4 ±0.3 ±0.6 ±0.4) x 10~ 4 , 
(1.8 ±0.3 ±0.3 ±0.2) x 10" 4 , (2) 



respectively. 

Since the masses of P*j are just below the pWff threshold and the corresponding widths 
are around few 10 KeV [ljj, both parity-even mesons could only decay through isospin 
violating channels to Dir and D*n. Due to this reason, it becomes the main problem how to 
explain the low masses and narrow widths for D* sJ states. Unlike D* sJ cases, however, there 
are no any suppressions on their decays to Dtt or D*tt although the measured masses of 
P** are slightly different from the predictions of theoretical models. It is believed that the 
properties of D** could be described by current theoretical models with some improvements. 
If so, based on the concept of the normal quark model, we could further understand the 
nature of D** in B decays. 

To handle the hadronic effects for B — > D**n decays, we use the factorization formalism, 
called perturbative QCD (PQCD) approach, which is based on factorization theorem and 



the transition matrix element is described by the convolution of hadron wave functions and 
the hard ta n d QQ. The wave tactions in principle can he extracted by expetaenta! 
data or determined by QCD sum rules or lattice calculations. The hard kernel is related 
to the hard gluon exchange and high energetic fermion propagator, which are all calculable 
perturbatively At the limit of heavy quark, in order to guarantee that color transparency 
mechanism [20] is satisfied, i.e., no soft gluon exchange between the final states, we need the 

it n 

hiera rchy of mg >> mo** >> A with A ~ m^-mj, ~ mo** — m c |2Jj. It has been shown by 
Ref. [22j that with the same QCD approach, the calculated results on P — > Dn decays are 
consistent with the current observations j^J. Consequently, one expects that PQCD could 
be also applied to the P-wave meson production in P decaying processes. By the study, we 
should know more properties related to P-wave mesons. 

The paper is organized as follows. We investigate the characteristic of D** and model their 
amplitude distributions in Sec. |Hj In terms of PQCD approach, we derive the factorization 
formulas for each B — > Dq7t decays in Sec. IIHI We present our results in Sec. IIVI Finally, 
we give a summary in Sec. El 

II. DECAY CONSTANTS AND WAVE FUNCTIONS OF D** 

In order to study the production of the scalar meson D** in B decays, we immediately 
have to face two questions. The first one is how to write down the hadronic structures and 
model the wave functions of D**, and the second is what the values of decay constants of 
D** are. In PQCD, since the wave functions belong to nonperturbative objects and are 
universal, we can directly apply the wave functions of B and 7r mesons, which have been 
discussed in other P-meson decays, such as P — > 7r£y, B — > tttt etc. Therefore, the hadronic 



structures of P and 7r mesons can be described by 



24. 



25] 



1 ^ 



(0|&((V(*),|fl,pi) = -== / dxe-^Wh + rriB^^Btx) 

V ^ V c JO 



(0|«(0)^(^|7r,p 3 > = -7^= / dxe-**»-*\\fa]Mx) 

\/21\ c J i 

+m° w [I} lj ^(x) + m<i[fi_fi + - J]^(*)] , (3) 

where n_ = (0, 1,0-r), n + = (1,0, 0y) and x is the momentum fraction of the light parton 
inside the corresponding meson. $ n and are the twist- 2 and twist-3 pion wave func- 
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tions, related to the distribution amplitude of nonlocal operator u n fs'y fl d(z) and associated 
with U75G?(z) (u^cr^d^z)), respectively. We note that is the so-called chiral symmetry 
breaking parameter and is equivalent to m\/(m u + m^). 

To determine the structures and distribution amplitudes of D**, we need to know their 
properties. For simplicity, we only concentrate the discussion on the scalar meson of Dq. 
The similar analysis can be applied to other charmed P-wave mesons. As usual, the decay 
constant of Dq is defined as 

(0\dc\D*,p 2 )=m D *J D * o . (4) 
By using the equation of motion, we obtain another identity 

(0|d7 M c|D3,p 2 ) = fD- P2^ (5) 

with fm = /ra(w c — Tnd)/ m D%, in which m c (d) are the current quark mass of c(cf)-quark. 
From the above equation, we see that if the considering case is light scalar meson, the 
corresponding transition matrix element will become small. This is the reason why the 
decay constant of the light scalar meson for vector current is small. In order to satisfy the 
conditions of Eqs. (jU and ©, the hadronic structure of Dq is adopted to be 

(D*Q,p 2 \d(0),c(z) l \0) = -^^ 1 «^-{[^] /j $ D1 (x)+m DS [/] / ,$ D2 (a;)}, (6) 



with the normalizations, 



The value of decay constant Jd* is the crucial part for concerning whether Dq production 
is interesting or not. To estimate the magnitude of fo*, we need the help with the scalar 
meson ^(1420), for which the decay constat has already been estimated in Ref. j^. As 
mentioned early, the scalar meson generally satisfies the identity 

(m qi - m q2 )(0\q 2 qi\S) = f s m 2 s , 

where m qi , ms and f$ are the current quark mass, the S'-meson mass and its decay constant 
of vector current, respectively. If we assume (0\d s|-K"q(1420)) ~ (0|<i c|Z?q), from the above 
equation, we can obtain f D * = f K * ■ m 2 K ,/m 2 D * ■ (m c — m d )/{m s — m d ). With the values 
of f K * ~ 34 MeV 3, m c = 1-5 GeV, m s = 150 and m d = 8.7 MeV, we get f D * w 130 



MeV. This value is close to the result in Ref. |27(, calculated by relativistic quark model. 
Finally, from Eq. (J5J) we have fo* = ^LD*/{f^c — Tn d ) ■ fo* ~ 200 MeV. It is known that Kq 
is composed of a two-quark state. Thus, it is interesting to have the similar decay constants 
between the scalar Dq and pseudoscalar D s . 

To obtain the shapes of Dq wave functions qualitatively, we need to employ the concept 
of the heavy quark limit. According to Eq. ©, we see that is the distribution 

amplitude of the nonlocal operator d'j fl c(z) while $02(2) is associated with dc{z). By 
the equation of motion, we straightforwardly find that the difference between §di(x) and 
&D2{x) is order of A/mn* ~ i m D* — mc)/ m D*- Hence, if we set m^* ~ m c , we can get the 
information of ~ &D2{%)- Furthermore, in order to satisfy the identities of decay 

constants defined by Eqs. (j3J) and (j5J), the simplest forms for both wave functions can be 
modelled by $£>j oc x(l — x) + a,ix(l — a;)(l — 2x) in which a, are free parameters. Since 
the second term is antisymmetric while x is replaced by 1 — x, we can easily conclude that 
this term will not change the normalization of the wave function. Therefore, we could use 
it to control the shapes of the wave function. It is worth to mention that since we consider 
Dq to be a P-wave state, the size of Dq is believed to be bigger than that of particle in the 
S'-wave state. In order to avoid that Dq becomes oversize such that the mechanism of color 
transparency is breakdown, like the 6-dependence on the wave function of the B meson, in 
which b is the conjugate variable of the parton transverse momentum, we also introduce the 
intrinsic 6-dependence on Dq. To satisfy Eqs. (J3J and (jSJ), the final simplest shapes of the 
wave functions are expressed as 



®D2(x,b) 



f. 



2y/2N c 



Id* 



6x(l — x) 



m c - m d 



+ a D *(l - 2x) 



2y/2N c 



6x(l -x)[l + b D *(l - 2x)] \ exp 



exp 



UJ 



2 h 2 
D*° 



(7) 



where u^*, o>d* and bp* are the unknown parameters. Although bp* is a free parameter, 



[m D * - 



it can be chosen such that the Dq meson wave function has the maximum at x 
m r ) /mn* ~ 0.35 for m c = 1.5 GeV. As to the value of an*, we refer to the case of ^(1410) 
2j|. By assuming that a D *m D */(m c - m d ) ~ a K *m K */(m s - m d ) fa 75/ f K *, the order of 
magnitude of am is estimated to be around 1.2. 
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III. FACTORIZATION FORMULAS 



Since the considered decays B — > D**ir correspond to the b — > cud transition, we describe 
the effective Hamiltonian as 

Hes = [CiifJ^daUpcpba + C 2 {[l)d a U a Cpbp\ , (8) 

where g^g^ = g«7 M (l — ^a)?/?, o;(/9) are the color indices, V c = V* d V c b is the pro duct of the 
CKM matrix elements 30], and Ci^{fj) are the Wilson coefficients (WCs) [3JJ. With the 
light-cone coordinate, the momenta of various mesons and the light valence quarks inside the 
corresponding mesons are assigned as: p\ = m B /\/2(l, 1, 0y), k\ = m B /\^2(xi, 0, kir)] p 2 = 
m B /V2(l, r 2 2 , T ), k 2 = m B /V2(x 2 , 0, k 2T )- p 3 = m B /V2(0, 1 - r 2 2 , T ), k 3 = m B /V2(0, (1 - 
rl)x 3 , k 3T ), with r 2 = m D * dj /m B . As usual, we use 

T = ^f^\V c \M? (9) 

to describe the decay rates of B — > D**tt, in which P c = \p 2z \ = \p 3 z\ ~ Tn B (l — r"^)/^ is the 
momentum of the outgoing meson, A is the decay amplitude and its value depends on QCD 
approaches. Since the hadronic structures of the tensor meson haven't been derived yet and 
so far they are not definite, we study the problem elsewhere. Although D[ and Di are the 
vector mesons and carry the spin degrees of freedom, only longitudinal polarization has the 
contribution since one of the final states is a pseudoscalar. Therefore, the deriving formulas 
for B — > DqTt are also proper to the final states with one vector and one pseudoscalar mesons. 
In this paper, we only concentrate on the production of Dq. 

In terms of the effective interactions, we see that different decaying processes involve 
different topologies. To be more clear, in the following we analyze each of B — > D^tt decays 
separately. 

1. — > Dq~ti + decay : 

There are two topologies in this decay, emission and annihilation diagrams. The former is 
color-allowed but the latter belongs to color-suppressed. The corresponding flavor diagrams 
are illustrated by Fig. [T] Hence, the decay amplitude of B d — > Dq~h + can be expressed by 

A{B d -> D* -n + ) = UF BD * + M BD * + f B F a + M a , (10) 

where f n and f B are the decay constants of 7r and B mesons and the related contributions 
are the factorizable emission and annihilation topologies, respectively. The remains denote 
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FIG. 1: The topologies (a) [(c)] factorizable emission [annihilation] and (b)[(d)] nonfactorizable 
effects for the decays B<i — > Dq~it + . 



nonfactorizable contributions. With factorization theorem and hadronic structures of Eqs. 
(JSJ) and the hard amplitudes {F} and {M} are formulated as 

/*1 /*oo 

Fbd* = 8nC F m 2 B / dxidx 2 \ b 1 db 1 b 2 db 2 &B(xi,b 1 ) 
Jo Jo 

x{$Dl(z2, &2)^D ( * (4 1) ) /i e(^l, Z 2 , &1, ^) 

+ (r c ^i(x 2 ) +2r 2 $ J , 2 (x 2 ))f| )S (t^ )h e {x 2 , x^hM)} (H) 

dx\dx 2 dx 3 I bidbib 3 db^B{xi,bi) 
Jo Jo 

X$ n (x 3 )$ m (x 2 , 6i){ - (x 2 + ^^(^/^({x}, h, b 3 ) 
^l-x 3 )8' 2 ,{tf)hf{{x}MM)}, (12) 

/•l POD 

F a = 87fC F m 2 B / dx 2 dx 3 / b 2 db 2 db 3 ^ D1 (x 2 ,b 2 )^ n (() 
Jo Jo 

{ ~ x 3 ^(t^)h a (x 2 ,x 3 r] 3 ,b 2 ,b 3 ) + x 2 £ 2 a (t {2) )h a (x 3 ,x 2 r] 3 ,b 3 ,b 2 )^ (13) 
M a = levrC^mlv^ / d[x] / [6] $ B , &i)$di(x 2 , & 2 )<MC) 

x{x 3 ^ 1 (4 1) )4 1) ({x},&l,&2)-^ 2 (4 2) )4 2) (W,&l,&2)}. (14) 

The hard functions h e ^,a,f), related to the propagators of exchange hard gluon and internal 
quark, are described by 

h e (xi,x 2 ,bi,b 2 ) = S t (x 2 )K (y/xiXirriBh) 

x [0(&i -b 2 )K {^/x^m B bi) I {^fx~ 2 m B b 2 ) 
+9(b 2 - bt)K {y/x~2~m B b2) h {\fx~2rn B bx)\ , 
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tq\x 1 ,x 2 ,x 3 ,b 1 ,b2) 



[9(bi - b 2 )K (yJxix 2 m B bi) I (y/x 1 x 2 m B b2) 
+9(b 2 - bx)K Q {^/xiX 2 m B b 2 ) I {■ s /x 1 x 2 m B bi)\ 
Ko{DjmBb 2 ) 



x 



for D) > 



ITT TT 



\D 2 \m B b 2 ) for D] < 



h a {x 2 ,x 3 ,b 2 ,b 3 ) 



StC^s) ^ J {\/x2X 3 m B b 2 ) 

9{b 2 - b 3 )H { Q ] ( v /a^m B 6 2 ) {\ / x 3 'm B b 3 
+6(b 3 - b 2 )H { Q ] (y^m B b 3 ) J {y/x^m B h) 



x 



^ ) ({x},& 1 ,6 2 ) 



7T 
I — 

2 



x 



6>(6i - & 2 )#d (^x 2 x 3 rj 3 m B b 1 ) J {^/x 2 x 3 r] 3 m B b 2 ) 
+6(b 2 - h)!!^ ( v /x 2 x 3 7] 3 m B b 2 ) J {^/x 2 x 3 r] 3 m B bi) 
KoiFjuiBh) for Ff > 



fif«( v /|F/|m B & 1 ) forF 2 <0 

with Df = xix 2 - x 2 x 3 7] 3} D\ = xix 2 -x 2 (l-x 3 )rj 3 , = (x 1 - x 2 )x 3 r] 3 , F 2 2 = x 1 +x 2 + (1 - 
Xi — £2)^3773, 773 = (1 — rf) and r c = m c /m B . The threshold resummation effect is expressed 
to be = 2 1+2c • T(3/2 + c)[x(l - x)) c /{^/ttT(1 + c)), with c « 0.35 [32]. The evolution 

factors £ l D *(£p*) and £ l a {£'a) are defined by 



Nr. 



a s (tf)exp[-S B - S D * - S n ] b , 2=bl , 



N n 



a s (tf)exp[-S B - S D * - S n } b , i=b2 , 



where the exponents Sm (M = B, Dq, 7r) are the Sudakov factors. From above equations, we 
see clearly that the emission contributions are color-allowed and dictated by effective cou- 
pling of C 2 + Ci/N c , while the annihilation contributions are color-suppressed and governed 
by C\ + C 2 /N c . t < f > da j denote the hard scales of the involving diagrams which are expected 
to be of 0(y/Am 2 B ) ~ 1.6 GeV in average and the criteria to determine them are adopted 



9 



to be 



max( v /xim B , 1/bi, l/6 2 ), = max(y/xim B , 1/h, l/b 2 ), 
max( v /xix 2 m s , ^Ljm B , l/b u l/& 3 ), 

max( v /i^m B , l/6 2 , l/6 3 ), t^ 2) = max( v /i^m B , l/6 2 , I/&3), 
max( 1 /i^m B , J F?m B , 1/h, l/b 2 ). 



(15) 



Since we deal with the hadronic effects of the B decay by considering six-quark simultane- 
ously, at lowest order in strong interaction, besides the renormalization group (RG) running 
from mw to m B scales in the /i-scale dependence of WCs, we still need to consider the run- 
ning from m B scale to the hard scale t^ da j. which indeed dictate the scale of the B meson 
decay. Hence, in our consideration, the hard scales for WCs are determined by Eq. ()15j) 
rather than at m B or m B /2 scale. In the formulations of Eqs. pijl — (|14 )1 . we have dropped 
the terms related to r\ (r c and r|) for the right-handed (left-handed) gluon exchange of Fig. 
[T] Compared to leading power, which isn't suppressed by l/m B , they all belong to higher 
power effects. 

2. B d Dg°n° decay : 




FIG. 2: The topologies (a) factorizable emission and (b) nonfactorizable effects for the decays 



B d D*°ir°. 



In this decay, the involving annihilation contributions are the same as the decay Bd — > 
Dq~h + but the emission topologies become color-suppressed, illustrated by Fig. EJ Due to 
the neutral pion meson being described by (uu — dd) / a/2, the sign of emission topologies is 
opposite to that of annihilation topologies. Therefore, the decay amplitude is written as 



A(B d 



1 

7! 



f D *F B7V -M B7T + f B F a + M a 



(16) 
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With the same approach and power counting for the Bd — > Dq 7r + decay, the relevant hard 
amplitudes F Bn (M Bn ) can be derived as 

/•I POO 

Fbit = 8irCFm B / dx\dx 3 / bidb 1 b 3 db 3 ^ B (xi,bi) 



x 



;i + x 3 )%(x 3 ) + r ff (l - 2x 3 ){<ffi{x 3 ) + ^(x 3 ))lf 3( f (3)) 



x/i e (xi,x 3 ?73,6i,63) + 2r 7r $£(x 3 )£ 7r (4 } )h e {x 3 ,xir] 3 ,b 3 ,b 

M Bn = l6irC F m B \J2N C / dx x dx 2 dx 3 l b 1 db 1 b 2 db 2 $ B (x 1 ,b 1 )§ D1 (x 2 ,b 2 

Jo Jo 

x | - {x 2 + x 3 )$ 7 r(x 3 ) + ?vr 3 ($£(x 3 ) + ^(ar 3 )) 
x£?(tf)h® (x 1 ,x 3 r ]3 ,x 2 ,b 1 ,b 2 ) 
+ [(1 - 12)^(2:3) " rvr 3 ($£(a; 3 ) + $£(x 3 )) 
x^ 4 (4 4) )/ i ? ) (xi,x 3% ,x 2 ,6 1 ,6 2 )}. 

The evolution factors E\{£'£) are defined by 



(17) 



(18) 



Ci(tf ) + 5^ija.(t?> )exp[-S B - 5.] 



C 2 (t 



a s (t { * ] ) exp[-S B - S D * - S n 



63=61 • 



From above equations, due to the appearance of C\ + C 2 /N c , we know that Bd 
is color-suppressed process. We note that although nonfactorizable effects are also color- 
suppressed, since C 2 /N c could be larger than C\ + C 2 /N c , the nonfactorizable effects play 

a important role in this kind of color-suppressed processes. In fact, the same thing also 

_ 

happens in the decay Bd — > D°7i° with D° being charmed pseudoscalar [22]. The hard 
scales are determined by 

tf ] = max(y/x 3 r) 3 m B , 1/h, l/b 3 ) , = max(y/xir] 3 m B , I/61, 1/63) , 

A (3) 



ty = max(^JxiX 3 r] 3 m B , \J D\m B , l/&i, l/6 2 ), 
4 4) = max( v /xix 3 77 3 m B , yD|m B , l/&i, l/6 2 ), 

with D| = (xi - x 2 )x 3 7] 3 and D| = (x x + x 2 )r| - (1 — xi — x 2 )x 3 i] 3 . 
3. B+ -»■ D*°7r+ deca?/ : 

In this decay, there are no annihilation contributions and new topologies involved. The 
corresponding flavor diagrams are the same as Fig. da) and (b) and Fig. |21 Hence, we can 



11 



immediately write the decay amplitude as 

A(B+ - D*%+) = f n F BD * + M BD . + f D *F Bn + M B ,. 
The hard amplitudes {F} and {M} are the same as Eqs. (ITT| . (ITIZJ) . (IT7I) and (fT8l. 



(19) 



IV. NUMERICAL ANALYSIS 

In our calculations, we adopt the B- meson wave function $^ to be 
§ B (x,b) = N B x 2 (l - x) 2 exp 





fxm B s 


) 2 ^ 2 j 


[-51 







(20) 

where N B can be determined by the normalization of the wave function at b = and uj b is 
the shape parameter. Since the 7r-meson wave functions have been derived in the framework 
of QCD sum rules, we display them up to twist-3 directly by |25[ 

3/. 



: x(i 



2V2iVe 



2^/2iv; 

with the Gegenbauer polynomials 



1 + 0.44C 2 3/2 (2x - 1) + 0.25C 4 3/2 (2a; - 1) 
1 + 0.43C 2 1/2 (2x - 1) + 0.09C 4 1/2 (2x - 1) 
1 - 2x) [l + 0.55(10x 2 - 10z + 1)] , 



c 2 3/2 (0 



Cl /2 (0 = ^(35£ 4 -30^ + 3), 



(5f -1), C 4 3/2 (0 = ^(21^ 4 -14f + 1 
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After the wave functions of 7r meson are determined, the unknown u B can be fixed by decays 
such as B — > mi. Consequently, the remaining uncertain parameters are the wave functions 
of the Dq meson. 

To obtain the numerical results, the values of theoretical inputs are chosen as: uo B = 0.4, 
f B = 0.19, U = 0.13, f D * = 0.20, m B = 5.28, m D * = 2.29, and m° = 1.4 GeV. With these 
values, we get the form factor F B ~' 7T (0) = 0.3. In addition, the values of B — > Dq form 
factor with some variances in a^* and ujd* are also shown in Table HI It is interesting that 
the form factor of B —>■ Dq decay is much smaller than that of B — ► D decay, which is 



calculated to be around 0.57 



221 ]. We also find that our results are a little bit larger than 

331 ] . we 



those calculated by ISGW2 model [291] . According to Wolfenstein's parametrization 
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TABLE I: The values of B — ► Dq form factor with bp* = 0.5 and some variances in ao* and wd*- 





a D * = 0.7 


a D * = 0.9 


a D * = 1.1 


0.5 


0.29 


0.30 


0.31 


0.6 


0.24 


0.25 


0.26 


0.7 


0.21 


0.22 


0.23 



take A = 0.82 and A = 0.22 for the CKM matrix element V c b = AX 2 . Hence, in terms of 
our deriving formulas and by fixing cld* = 0.9, bo* = 0.5 and oj£>* = 0.6, the magnitudes of 
the hard amplitudes are shown in Table |HJ From the table, we can see clearly that except 

TABLE II: The values of hard amplitudes (in units of 10 -3 ) with fixing ap* = 0.9, bp* = 0.5 and 
u D * = 0.6. 



fnFBD* 


M BD * 


f B F a 


M a 


Id*Fb-k 


M B n 


36.4 


10- 2 (1.0-i3.2) 


-0.06 - i0.08 


-1.85 - z3.14 


-7.11 


7.56 - il0.95 



Mbd^, the nonfactorizable effects of color-suppressed process are comparable to factorizable 
contributions; even in annihilation topologies, the contributions of the former are much 
larger than those of the latter. With fixing bp* = 0.5 and ujd* = 0.6 GeV and taking some 
different values of od* , the decay BRs of B — > DqTt are displayed in Table IIHI We also show 
the BRs with fixing an* = 0.9 and bp* = 0.5 and some variant values of u)£>*. From both 
tables, we know that with proper values of parameters, the calculated BR of B + — > Dq°7t + 
is consistent with the BELLE's observation. It is worth to note that the predicted BR 
of Bd — > Dq°ii is one order of magnitude smaller than others. The phenomenon can be 
understood by noticing that, as shown in Table ILT1 the value of Fb-k is very close and opposite 
in sign to the real part of Mb-k such that there is a strong cancellation in Eq. (|T6*j) . As a 
result, we get the small BR in the decay B^ — > D^n . That is, the annihilation effects are 
significant in B — » D^n . 

As stated before, although we only study the decays B — > -DqVt, we still can estimate the 
BRs of B — > D^tc. Since only the longitudinal polarization has the contributions, except 
the decay constants, we expect that the involving wave functions of D± should be similar to 
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TABLE III: The BRs (in units of 10 4 ) with fixing br>* = 0.5, ujd* = 0.6 GeV and various values 
of an* ■ 



a-Ds 


B + -» D^°TT + 


B d - L»*-7r+ 


B d - ^o* ^ 


1.1 


9.75 


8.25 


0.17 


0.9 


9.34 


7.68 


0.19 


0.7 


8.98 


7.13 


0.21 



TABLE IV: The BRs (in units of 10 4 ) with fixing cld* = 0.9 and bp* = 0.5 and various values of 





B+ -> L^ 7T+ 




5 d -> 


0.5 


13.79 


10.7 


0.26 


0.6 


9.34 


7.68 


0.19 


0.7 


6.28 


5.55 


0.15 



Dq. By neglecting the difference in phase space, the BRs of B — > D-pw could be estimated 
by BR(B -> D { ^)/BR(B -> D*tt) ~ (f^/fn*) 2 . If / D w « /d S , the BRs for producing 
axial vector mesons Dj are close to that for the scalar Dq. The tendency is consistent with 
BELLE's observations, shown in Eq. (J2J). 

V. SUMMARY 

We have studied the properties of P-wave mesons in B decays in terms of Dq. By taking 
the concept of the heavy quark limit, we have obtained some information on the shapes 
of Dq wave functions. According to the wave function of _K"q(1420), we can determine the 
proper value for the parameter ap* in By the physical argument, the unknown 

parameter bo* can be chosen so that the maximum of §D2{x) locates at x ~ 0.35. We have 
found that with a suitable value of uj£>*, our result on BR(B + — > Dq°tt + ) can fit BELLE's 
measurements. Hence, the calculated BRs for Bd — > Dq~tt + and Bd — > Dq 7t° decays can 
be viewed as our predictions. Finally, if we regard that the longitudinal wave functions of 



14 



are the same as £#° and assume that f D (,) ~ /d*, we expect that the differences of 
BRs among them are not significant. The more accurate predictions rely on more definite 
values of decay constants as well as other unknown parameters. 
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